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ABSTRACT 
> 

Cf^ , Using multi-frequency spectra from TeV blazars in quiescent states, we obtain 

^ [ the physical parameters of the emission region of blazars within the framework 

O ! of the one-zone synchrotron self-Compton (SSC) model. We numerically 

calculate the steady-state energy spectra of electrons by self-consistently 
r-| ! taking into account the effects of radiative cooling with a proper account of 

5^1 the Klein-Nishina effects. Here electrons are assumed to be injected with a 

O I power-law spectrum and to escape on a finite time scale, which naturally leads to 

^ [ the existence of a break energy scale. Although we do not use time variabilities 

^ ' but utilize a model of electron escape to constrain the size of the emission 

region, the resultant size turns out to be similar to that obtained based on time 
r> I variabilities. Through detailed comparison of the predicted emission spectra 

c3 : with observations, we find that for Mrk 421, Mrk 501, and PKS 2155-304, the 

energy density of relativistic electrons is about an order of magnitude larger 
than that of magnetic fields with an uncertainty within a factor of a few. 



Subject headings: BL Lacertae objects: general-gamma rays: theory-radiation 
mechanisms: nonthermal 



1. INTRODUCTION 

Blazars comprising of BL Lac objects and optically violent variable quasars are 
characterized by rapid time variation of the energy flux, large and variable polarization. 
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and featureless continuum spectra (e.g., Urry & Padovani 1995). These characteristics 
are considered to be the result of beamed emission from relativistic jets seen end-on (e.g., 
Blandford & Konigl 1979). The discovery of strong 7-ray emission from blazars in the 
GeV band by EGRET on Compton Gamma Ray Observatory is one of the most important 
issues for active galactic nuclei (AGNs), because more than 60 AGNs detected by EGRET 
are all blazar type and no identifications as other types of AGNs such as Seyfert galaxies 
have been reported (Mukherjee et al. 1997) except for a probable detection from the radio 
galaxy Centaurus A (Hartman et al. 1999). Multi-frequency observations have revealed 
that broad band continuum spectra of blazars consist of two components; the low energy 
component from radio to optical/UV band sometimes extending to X-ray band is by 
synchrotron radiation, while the high energy component from X-ray to 7-ray is due to the 
inverse Compton scattering of soft photons (e.g., Kubo et al. 1998): various soft photon 
sources have been proposed ranging from synchrotron photons to photons from accretion 
disks, either direct or reprocessed (e.g., Sikora, Begelman & Rees 1994; Inoue & Takahara 
1996; Dermer & Schlickeiser 1993; Blandford & Levinson 1995; Ghisellini & Madau 1996). 
One of the most important aspects of multi-frequency observations of blazars is to probe 
the energetics of relativistic jets. From this point, a few blazars from which TeV 7-rays 
have been detected (Mrk 421, Mrk 501, PKS 2155-304, and lES 2344+514) are especially 
important, because these TeV blazars are relatively less luminous and pure synchrotron 
self-Compton model can be best applied. 

Up to now, source parameters of TeV blazars have been estimated in a variety of ways 
(e.g., Bednarek & Protheroe 1997; Tavecchio et al. 1998; Kataoka et al. 2000). With regard 
to the basic energetics, however, surprisingly little attention has been paid to the energetics 
of electrons, which should be the central concern for the theoretical understanding of 
the production and bulk acceleration of relativistic jets. Based on an analytic estimate 
from TeV blazar observations, one of the present authors argued that relativistic electrons 
dominate over magnetic fields in energy densities in relativistic jets of blazars (Takahara 
1997). The purpose of this study is to estimate the energy densities more quantitatively 
using a numerical code which self-consistently solves for electron spectra suffering from 
injection, escape, and radiative cooling and photon spectra with proper account of the 
Klein-Nishina effects (e.g., Mastichiadis & Kirk 1997; Li & Kusunose 2000). Although 
recent observations show that even in low activity states 7-ray spectra extend above 1 
TeV (e.g., Aharonian et al. 2001), in this paper, we restrict our attention to the quiescent 
states and treat the 7-rays below ITeV in the first step neglecting the correction for the 
absorption of TeV 7-rays due to Cosmic Infra- Red Background (GIB). We shortly discuss 
GIB absorption effects in §|^. In a future research, we will separately examine this issue 
including the flaring states where 7-ray spectrum clearly extends up to 10 TeV. 
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The rest of the paper is organized as follows. In §|^, we introduce the synchrotron 
self-Compton (SSC) model, and we describe our numerical treatment for solving the kinetic 
equations of photons and electrons. We also show the relation between model parameters 
and typical observables to help the search for the parameter set of the best fit model 
in numerical calculations. In §^, we discuss the analytic estimation of the ratio of the 
energy density of relativistic electrons to that of magnetic fields. In §^, we show the 
numerical results of spectral fitting applied to three TeV blazars, i.e., Mrk 421, Mrk 501, 
and PKS2155-304. (1ES2344+514 is omitted because less data are available at present.) 
Finally in §^ we summarize our main results and discuss some related issues. 



2. ONE-ZONE SSC MODEL 

2.1. Basic Assumptions 

Non-thermal emission from TeV blazars is divided into two components, i.e., low 
energy synchrotron component extending from radio to X-rays, and high energy inverse 
Compton component extending from hard X-rays to TeV 7-rays. Here the seed photons of 
inverse Compton scattering are the synchrotron photons in the same emission region. This 
SSC model has been very successful in describing the observed multi-frequency spectra 
(e.g., Jones et al. 1974; Maraschi, Ghisellini, & Celotti 1992). 

Further assumptions used in the present work are that (1) emission region is one zone 
with a characteristic size R and is moving at a relativistic speed (3 in units of the light 
speed and that (2) both relativistic electrons and photons are isotropic in the source frame. 
Beaming (Doppler) factor is given hj 6 = l/[r(l — jScosO)], where 6 is the angle between 
the line of sight and the direction of the relativistic jet and F is the bulk Lorentz factor 
of the emission region in the jet. When the observer lies within the angle of ^ ~ 1/F, we 
obtain F ~ 5. The Hubble constant is assumed to be 75 km s~^ Mpc~^. Throughout this 
paper we use these approximations. 



2.2. Numerical Approach 

Most of previous calculations are either semi-analytic or done without the inverse 
Compton process self-consistently. In our numerical code, to obtain the consistent spectra 
of photons and relativistic electrons, we calculate the kinetic equations of electrons and 
photons self-consistently including the exact inverse Compton process within the continuous 
energy loss approximation. 
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The kinetic equation describing the time evolution of the electron distribution is given 
by 

Q-^ + — = -Q- [(7syn + 7ssc)ne(7, ^j] + Qe,mj(7, t), (1) 

where 7 is the electron Lorentz factor and 11^ is the electron number density per 7; 7gyn 
and 7ssc are the cooling rates of synchrotron and inverse Compton emission, respectively; 
te,esc is the effective escape time of the electrons, which is identified as the time scale of the 
adiabatic expansion loss (Mastichiadis & Kirk 1997). For simplicity, in all the following 
numerical calculations, we set te,esc = 3tdyn, where tdyn = R/c. Next we adopt an injection 
spectrum 

Qcinj = ge7"'e-^/^-^'' for7^i„<7, (2) 

where 7max and 7min are, respectively, the maximum and minimum Lorentz factors of the 
electrons, qe is the normalization factor, and s is the power-law index. As for the injection 
mechanism, we implicitly assume the first order Fermi acceleration (e.g., Blandford & 
Eichler 1987). The synchrotron emissivity and absorption coefficient are calculated based 
Robinson fc Melrose (1984]| and [Crusius fc Schlickeiser (1986)| . In the calculation of the 



on 



inverse Compton scattering, we use the exact Klein-Nishina cross section and scattering 
probability of [Jones (1968)| and [Coppi fc Blandford (1990) 



The kinetic equation of photons is given by 
(9raph(e,t) nph(e,t) 



dt t 



ph,esc 



nic(e,t) +nsyn(e,i) (3) 



where e is the dimensionless photon energy normalized by r/igC^ with m^ being the electron 
mass, Uph is the photon number density per unit energy e, tph,esc is the escape time of 
photons from the emission region, which is taken as tph,esc = R/c in the optically thin limit. 
hic{e,t) and nsyn(e, i) are the production rate of inverse Compton and synchrotron photons 
per unit energy e, respectively. 

In order to obtain quiescent state spectra, calculations are done up to IStdyn, which is 
long enough to reach a steady state. The physical quantities in the source frame can be 
converted to those in the observer frame using the relations such as Co = es(5/(l + z) and 
dto = dts{l + z)/6, where subscripts o and s express the quantity in the observer and source 
frame, respectively, and z is the redshift of the source. 

In this model, there are seven parameters to be determined by the comparison of 
predicted and observed photon spectra. They are: R the size of the emission region, B the 
magnetic field strength, 6 the beaming factor, 7max the maximum Lorentz factor, 7min the 



-5- 



minimuin Lorentz factor, q^ the injection rate of electrons, and s the power-law index of the 
injected electron spectrum. Among them, 7jnin is not easily constrained by spectral fitting 
and it is taken to be 10 in all the numerical calculations. Although 7inin hardly affects the 
radiation spectra, it is important for probing the energy and number densities of relativistic 
electrons and thus the matter content of the relativistic jets, and it has been a matter of 
debate (e.g., Reynolds et al. 1996; Wardle et al. 1998; Hirotani et al. 1999). We discuss the 
effect of changing the value of 7min in §|. 



2.3. Analytic Estimate 

Before we present numerical results, we describe some analytic estimates which provide 
useful insight into the physics behind the relationship between the model parameters and 
typical observables and we later examine quantitatively to what extent simple analytic 
methods are accurate. Using these relations, we constitute an analytic estimate of model 
parameters which are then used as a starting set of parameters for numerical calculations. 



2.3.1. Relation between the Model Parameters and Observables 

Among seven parameters, 7min is least constrained and taken to be 10 as was mentioned 
above. The index s is determined by the spectral shape of the synchrotron radiation at 
low energies; specifically, the energy index a between the radio and IR band is used to 
determine s by 

s = 2a + l. (4) 

Other five parameters, R, B, S, 7max, and qe remain to be determined. Basically, the 
luminosities and typical frequencies of synchrotron and inverse- Compt on components give 
four constraints. Remaining one can be taken to be the break frequency of the synchrotron 
radiation which corresponds to the break Lorentz factor of electrons, 7br, resulting from 
radiative cooling before escape. In principle, the break feature may appear in the Compton 
component, too. However, the spectral resolution of the present 7-ray observations is not 
good enough. Moreover, the Klein-Nishina effect makes the situation complicated. Thus we 
do not use the break frequency of the Compton component in this paper. 

To sum up, the five typical observables in the observer frame are: t'syn.o.max the 
maximum synchrotron frequency, t'syn.o.br the synchrotron break frequency, z/ssc,o,max the 
maximum frequency of the SSC component, -Lsyn,o total synchrotron luminosity, and -Lssc,o 
total SSC luminosity. Schematic pictures of a multi-frequency radiation spectrum and a 
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relativistic electron energy distribution are shown in Figures ^ and |^, respectively. The 
approximate solution of the electron kinetic equation (|1]) is 



Ml) = 


Qe'^e,escl 


for 


7min < 7 < 7br 


Ml) = 


~ ?e^e,escTbrT 


for 


7br < 7 < 7max 


provided that 7min < 7br- 









(5) 



Using the standard formula about radiation (e.g., Rybicki & Lightman 1979), we can 
obtain five relations between the model parameters and observables. Observed synchrotron 
frequencies from a single electron of Lorentz factor 7br and 7max are respectively given by, 

r 

Z/syn,o,br = 1.2 X IQ^^B^l—- (6) 

i. ~j~ z 

and 

i^syn,o,max = 1-2 X IQ^ B^^^^^-^^. (7) 

The maximum value of the observed SSC energy in the Klein-Nishina regime is 

''''^ssc,o,max '-^l7max''''c'-- 



1 + z' 

where h is the Planck constant and Ci < 1 is a constant representing the uncertainty of the 
Klein-Nishina effect which typically taken to be 1/3 here. It is to be noted that hvssc,o,ma.x 
is limited by the Klein-Nishina effect unless the beaming factor is extremely large (typically 
larger than about 100). This is understood as follows. The detection of TeV photons means 
that 7niax is at least greater than 10^'^/(5. If (5eseed,s7max ~ 10^^ ^V and 7niax > 10^'^/(5 are 
satisfied, we obtain eseed,o < 0.15^ eV, where egeed.s and esced.o are seed photon energy in 
the source frame and the observer frame, respectively. This means that the observed seed 
photon frequency is lower than X-ray band unless 5 > 100. 

In the case of s < 5/2, bolometric synchrotron luminosity in the observer frame is 
given by 

A-TT /? /"7max A. 

^syn,o = 47rDL^syn,o = ^^— ^^ / ^ ^TC7^MBge4,esc7br7~'~^C?7, (9) 

where mb is the energy density of magnetic fields, Dl is the luminosity distance, and -Fsyn,o 
is the total flux of synchrotron radiation in the observer frame. Note that in the case of 
s > 5/2, the luminosity from electrons with imm < 7 < 7br is larger than that from electrons 
with 7br < 7 < 7max and equation (^ is not valid. The bolometric Compton luminosity 
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is written in a similar way, by replacing ub by the energy density of soft photons and 
multiplying a suppression factor of C2 mentioned below. The energy density of synchrotron 
photons is given by 

-^syn.o = 7^ C^ Wsyn , (lU) 

when we set the photon escape time as R/c. It is important to note that because of the 
Klein-Nishina suppression only photons with energy less than meC^/7 in the source frame 
contribute to SSC luminosity. Here, we simply denote this suppression factor by C2. The 
ratio of the synchrotron and SSC luminosities is then given by 



-^syn,o ^B 



The break Lorentz factor is determined by the condition that te,esc is equal to tcoob 

R SnieC 

C^C 4(mb + C'2Msyn)0^T7br ' 



;ii) 



(12) 



where we set Cste^esc = ^dyn and assume C3 = 1/3 as has been noted in the previous 
subsection. This is a different approach from most of the previous work where the time 
variability constraint R ~ 5ct^ax is used. We think it more appropriate to avoid time 
variability constraint in the case of quiescent states, because the shortest time variability 
such as 15 minutes TeV flare (Gaidos et al. 1996) might be correlated to local regions such 
as a shock front (Kirk et al. 1998). 



2.3.2. Analytic Estimate of Physical Parameters 

In this subsection, we analytically estimate the model parameters using the typical 
observables of TeV blazars. Observed total flux, typical frequencies, and luminosity distance 
are scaled as 

(. -'^syiijO (. -'^ssCjO /, qn 

Jsyn -, ri — in ^ IT ' ^ssc -, n — in To ZT i K^'-'l 

lU ^^erg cm "^ s ^ 10 ^^erg cm ^ s ^ 



'^^syn,o,br '^syn,o,max '^ssc,o,max 



and 



''syn,o,Dr ''syn,o,max "ssCjOjinax ( ^ A\ 

''"'^"10^^' ''"^"^^ 1017Hz ' ''^^^^ 1026Hz ^ ^ 



d=-^^, (15) 



100 Mpc ' 



respectively. The numerical factors Ci, C2, and C3 normalized by 1/3 are denoted by Ci, C2, 
and C3, respectively. 

In order to express the model parameters in terms of the observables, first we solve 
five algebraic equations (|), (0), (§), (|ll]), and (|12D, to obtain five quantities R, B, 6, 7br, 



and 7max5 and then by inserting them into Eq. (||), we obtain q^. As a result, we obtain the 
typical values of parameters for TeV blazars as following: 

5 = 0.13 4/,Vs;c^^^/V^/V^/4.-.^/^rf^/^(l + zr crepes'/' G , (17) 

7_ = 2.7 X lOV-y Vil^r,-^/^-.-//V-^/V3/^rf-V2(l + z)cf/%'/^/' , (19) 
/? = 9.0 X lO^Vs^y Vi^^r^-^/V-//V-I/VJ/,^ci-V2(i + z)cf/W/'cr cm , (20) 



and 



QeT^l-'-'d^ = 0.015 /VV^^r,^/V,f.i^/>-V^d^/^(l + .)-3 

"'7br 

xcY/'c^'/%'/' cm--' s-\ (21) 



where 77 = (1 + Fssc,o/-f^syn,o)/2. Here we emphasize that these expressions give a complete 
set of model parameters in terms of the observables and that typical numerical values turn 
out to be similar to those obtained in various other ways (e.g., Tavecchio et al. 1998; 
Kataoka et al. 2000). In particular, the obtained size seems to be similar to that based 
on time variability, which means that the size of the emission region in quiescent states is 
compatible with that obtained by time variability constraint. It is seen that R and qe have 
a strong dependence on Umax and z/ggc, and that a variation of z/ggc by a factor of 2 leads to a 
variation of R, B, and qe by a few orders of magnitude. Despite this, quite robust estimate 
is possible for some of the quantities such as the ratio of the energy density of electrons to 
that of magnetic fields, as will be shown in the next section. 

Although our method is model dependent in that we use a simple model of electron 
injection, escape, and cooling, this model is quite general and has an advantage of 
self-consistent treatment of the break Lorentz factor of electrons. If one does not use this 
relation and tries to proceed based on observables alone, one needs to introduce the electron 
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number density instead of the injection parameter q^ and to use time variability constraint. 
Most of the previous work adopted such methods and searched for suitable parameters 
in a two dimensional parameters such as (logS, log 5) plane (e.g., Tavecchio et al. 1998; 
Bednarek & Protheroe 1997; Kataoka et al. 1999) by allowing some degree of uncertainties 
of time variability constraint. 

Next we check an additional constraint on the transparency of 7-rays against the 
intrinsic absorption, i.e., the optical depth of a 7-ray photon for pair production should be 
smaller than unity. For a 7-ray photon of the observed energy e^^o, the observed energy of 
the target photon et,o is about 

* ... = ^^(v^f. (22) 



1 + z e^,o 1 + z 

We approximate the number density of target photons by 






Then, the optical depth is given by (von Montigny et al. 1995): 



15 1 + z e^^pLt^oCTT 



9-9 X ^^-'^7J::^^^r f TT^r' (4)-^(l + -)• (24) 



lO^^erg s-i TeV V W^cuiJ ' 10 
If we use the typical observables, the optical depth is expressed as 

r,,(e,) = 9.8 x 10"^ ^f-lUrW^UlA^ + zr'c,c:,'c, ^^'° . (25) 

Thus we find that intrinsic absorption of 0.35 TeV photons is not so large for TeV blazars 
in quiescent states but it may be important for some cases. In [Q, we will show the value 
for numerical solutions in Tables 1, 2, and 3. 



3. Energy Density of Relativistic Electrons and Magnetic Field 



In this section, following the methods of [lakahara (1997)] , we show that we can 



analytically estimate the ratio of u-b to Ue fairly robustly. The energy density of relativistic 
electrons is given by 

/■7max 

tie = / im^c ne(7)ci7 

■'7miii 
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2 n^' - 4-1 



mm 



1.1 X lO-'f^^n'^^^'^^^^^dil + z)-\\cr'-'cf X g:^^t!. erg cm-3(26) 



-2 f 1/2^2,^1/2,^5/2 ,^-5^n ^ .^"2^5 ^7^2 ^-2 ^ ^Tmi. 



in the framework of the previous section. The energy density of magnetic fields is also given 
by 

«B = 6.3xlO-%n/s;c'Vbf^^il^s;crf(l + ^)"'c?cfc3-' crg cm-^ (27) 

Although both energy densities have fairly strong dependence on some observables, the 
ratio turns out to be as simple as 

— = 9 C2-'c,-' fl + ^^ ^ !r"^"Tit^ • (28) 

Note that the strong dependence on quantities such as z/gsc, z^max, and ci in Me and mb are 
canceled out. In addition, we can derive equation (^8]) in a more transparent way as follows. 
The total radiation power is simply given by 






max 



i^syn.o + ^ssc,o = -TiR^d'^m^c^q^ / 7-^+^7. (29) 

J "'7br 



Combining equations (0), (|Tl]), (pBf), and (^) we can easily obtain equation ( pS]) (The 
factor 9 in Eq.(|28|) comes from the normalization of C2 and C3.). This analytic estimation 
is quite useful to understand the relation between the typical observables of TeV blazars 
and the ratio of Uc/u-q. Needless to say, the ratio of 7br to 7max is obtained from the 
observed ratio of z/syn,o,br to t'syn.o.max (Eqs. H and 0). The ratio of 7br to 7min depends 
on the adopted value of 7min- When the synchrotron luminosity dominates over the SSC 
luminosity, -Usyn ^ mr and Uc/u-q oc /ssc//syn ^ 1- When the SSC luminosity dominates 
over the synchrotron luminosity, -Usyn ^ mb and Uc/ub oc fssc/fsyn ^ 1- Hence, as the value 
of /ssc//syn iucreases, the value of Us/ub increases. Very roughly, the equi-partition between 
electrons and magnetic fields corresponds to sources for the SSC luminosity equal to the 
synchrotron luminosity, if we ignore other numerical factors such as C2 and C3. According to 
equation (^), more realistic estimate indicates that when the SSC luminosity is equal to 
the synchrotron one, Ue/uB takes a value of an order of 10. 

Since the above estimate depends on several numerical factors, we should examine 
carefully these subtleties. As for the value of C2, the energy of seed photons for 7-rays 
of e^^o TeV is lower than about 10e~Q eV for a typical value oi 6 = 10. Because of the 
Klein-Nishina effects, synchrotron photons whose energy is higher than this value is not 
available to inverse Compton scattering. Thus C2 is expected to be less than 0.3. As for 
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e's, since electron escape may be identified with the expansion velocity in the downstream 
region of the shock, 1/3 is also a reasonable guess; higher values approaching 1 requires as 
rapid escape of electrons as that of photons, which seems to be difficult to realize. The 
typical ratio of 7max to 7br is 10. Even if this ratio is 100, u^ju^ decreases only by a factor 
of a few when s = 1. An increase in the value of 7min by a factor of 10 will decrease u^/u-q 
by a factor of a few for s = 2. Thus, only when both 7inin and 7max/7br are larger than 100, 
Mb can be comparable to u^- We should note that since 7min is expected to be comparable 
with the value of the bulk Lorentz factor in the shock acceleration picture, 7min = 10 is 
a reasonable choice when the jet material consists of mainly electron-positron pairs or 
when electrons and ions are separately thermalized for normal plasmas (see §^, for further 
discussion) . 

Thus, the only major uncertainty in estimating us/ue is the electron index s. In a 
previous work of P?akahara (1997) , he examined the case of s = 2, because this is the 



universal index expected for non-relativistic case (e.g., Blandford & Eichler 1987). In the 
present work, we allow this value to be less than 2, which better fits the emission spectra. 

When s is smaller, Uq/ub becomes smaller. In Figure ^ we show the resultant ratio 
Ue/uB as a function of 7br for several values of s. For simplicity, other parameters are fixed 
at /syn = 1, /ssc = 1, 7max = 1 X 10^ 7min = 10, Ci = 1, C2 = 1, and C3 = 1. Figure | shows 
that even when s is as small as 1.4, Ug is larger than ub, unless 7br is smaller than 10'^. 
Thus, the conclusion that u^ is about one order of magnitude larger than ub is fairly robust. 

Before we present detailed numerical results, we describe straightforward relations 

between the energy densities and power carried by relativistic jets. It should be noted that 

the energy density of i-th radiation component (the suffix i corresponds to synchrotron or 

SSC), Ui is related to the power Lj as 

4 
L, = -TiR^cUiV'^ , (30) 

where the factor 4/3 accounts for pressure of relativistic matter and we assume F ^ 1. It is 
worth noting that the observed luminosity is given by 

Li,o = ^TiR^cu.d'' (31) 

in the observer frame. Hence the actual radiation power is smaller than the observed 
luminosity by a factor of 5'^ assuming F ~ 5. Similarly, the Poynting power is given by 

Lpoy = ^ttR^cubT^ (32) 

and the kinetic power of relativistic electrons is 

4 

^e,km = -7ri?^CUeF^ . (33) 
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Thus, the discussion on the ratio u^/u-b is straightforwardly translated into the ratio 
Lc,km/Lpoy It should be stressed that -Le,kin takes account of only relativistic electrons. In 
other words, the contribution from thermal electrons which should constitute a reservoir for 
acceleration and that of protons, either relativistic or cold, are completely neglected. We 
emphasize that, if we take into account of these components, it is clear that kinetic power 
is more dominant at least in the emission region of TeV blazars. 



4. NUMERICAL RESULTS 

Here, we present numerical results by solving the kinetic equations and searching for 
best fitted model parameters for three TeV blazars (Mrk 421, Mrk 501, and PKS 2155-304). 
The search is made around the parameters determined analytically in the previous section. 
The observed values are taken from those compiled by [Kataoka (2000)| . In his paper, the 



multi-frequency observed spectrum is fitted with a polynomial function of the form of 
logi^uL^) = a + blogu + c(logz/)^ + dilogu)^, where a, b, c, and d are the fitting constants 
(Comastri, Molendi, & Ghisellini 1995). [Kataoka ( 2000 )| identified z/syn,o,br and t'syn.o.max as 



the frequencies where synchrotron luminosity reaches a half of its peak value; the lower 
one is identified as i^syn,o,br and the higher one is identified as t'syn.o.max- In the same way, 
z/ssc,o,max cau bc defined. However, since the spectra of high energy 7-rays are not so well 
covered, we do not put a heavy weight on his determination of z/ssc,o,max, but we simply 
assume t'ssc.o.max = 1 X 10^^ Hz as a starting point. It should be noted that these observables 
are not derived from the spectral fitting described in the previous section, but from the 
polynomial fitting, so that the chosen values of the model parameters above are expected 
to deviate from our definition by some factor. 

Moreover, taking the sparseness and uncertainties of the observed data into 
consideration, it is natural to think some range of uncertainties are included in the model 
parameters, too. Such uncertainties are also investigated in this section. Here, we explain 
the procedure of searching these parameter sets. The range of uncertainties is different for 
each model parameter. Among seven model parameters, the index s can be regarded as 
fixed. In contrast, R and qc are expected to have a large uncertainty, while 6 is relatively 
stable. The search for model parameters is made as follows. We first choose a certain value 
of ge- Then, we adopt a suitable value of R so as to reproduce the low energy spectrum of 
the synchrotron component, whereby some discrepancies will appear at the high energy part 
of the synchrotron component. Noting that the observed ratio of Lssc,o/-Z^syn,o is determined 
by the combination of B^5^R^, we can adjust the break feature of the synchrotron spectrum 
by adopting a suitable combination of B and 5. A slight adjustment of the value of R 
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is also made. Finally, the high energy end of the synchrotron component is adjusted by 
adopting a suitable 7max- In this step, slight adjustment of the model parameters which 
are determined in the previous steps is also made. The resultant spectral shape of the 
SSC component should be compared with observation, while the calculated SSC luminosity 
should match the observation. This comparison determines whether the chosen parameter 
set is allowable or not. The whole step is repeated starting from a different value of q^ 
to find the best fit model parameters and uncertainties. Although this is not a complete 
survey of the parameter space, we believe that this provides a reasonable estimate of the 
uncertainties because the most uncertain parameters are q^ and R. 



4.1. Best Fit Parameters 

4.1.1. Mrk 421 

Mrk 421 {z = 0.031) is a well known BL Lac object and the first identified source of 
TeV gamma-ray emission by Whipple Cherenkov telescope (Punch et al. 1992). 



Following [Kataoka (2000)| , the observables of Mrk 421 are chosen as follows; 



a = 0.3, I/syn,o,br = 2.5 X 10^^ Hz, I/syn,o,max = 1-6 X lO^^^ Hz, I/ssc,o,max = 1-0 X 10^^ Hz, 

-^syn.o = 8.6 X 10"^° erg cm~^ s~\ and -Fssc,o = 3.4 x 10"^° erg cm"^ s"^ Substituting these 
observables into the analytic estimate described in §2, we obtain 6 = 14.3, R = 7.9 x 10^^ 
cm, B = 0.31 G, 7br = 0.69 x 10^ 7inax = 1-8 x 10^ ge = 3.2 x 10"^ cm'^ s-\ and s = 1.6. 
Dotted line in Figure ^ shows the predicted spectrum obtained by using these analytic 
values. 

This first set of parameters in fact produces much higher luminosities than observations, 
which is not surprising as explained above. Then we try to search for the best fit parameters 
by changing model parameters to match the predictions with observations. Finally, we 
find a satisfactory set of model parameters as S = 12, R = 2.8 x 10^^ cm, B = 0.12 G, 
7max = 1-5 X 10^, ge = 9.6 X 10~^ cm~^ s~^, and s = 1.6. In this fitting procedure, we fix the 
reference value of s, because s can be determined with little error from the spectral shape of 
the low energy synchrotron emission. The thick solid curve in Figure ^ shows the predicted 
spectrum of Mrk 421 calculated from equations (^ and (^ for the best parameter set given 
above. The corresponding electron energy spectrum is shown by the thick solid curve in 
Figure 5. The ratio of Uq/u-q derived using this self-consistent numerical result is Uc/u-q = 5. 

The thin solid curve and the dashed curve in Figures ^ and ^ are for different sets 
of model parameters to examine the range of uncertainties in the values of the model 
parameters. The dashed curve is for the case that qe is 50 times larger than the best fitted 
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value. The thin sohd curve is for the case that q^. is 50 times smaller than the best fit 
value. These two sets of parameters may be regarded as marginally allowed, although the 
discrepancy in the TeV range is fairly large. The resultant values of the model parameters 
are tabulated in Table 1. As is seen, for the larger value of ge, the values of R and 7max 
become smaller, while the values of B and 5 become larger. For a smaller value of ge, 
the reverse is true. Table 1 shows the range of uncertainties for model parameters too. 
In particular, the ratio of u^/ub is uncertain by a factor of a few and the dominance of 
electrons over magnetic field is not changed. The values of 6, 7max, and B are uncertain 
by factors of about 2, 5, and 10, respectively. The size of the emission region is least 
constrained with an uncertainty by a factor of 30, but covers a reasonable range. These 
results can be roughly understood by combining equations (||), (|^, (||), and (|lT]). 



4.1.2. Mrk 501 

BL Lac object Mrk 501 {z = 0.034) is also a well known source of TeV 7-rays detected 
by Whipple Cherenkov telescope (Quinn et al. 1996). Observables of Mrk 501 are chosen 
as follows (Kataoka 2000); a = 0.4, i^syn,o,br = 6.3 x 10^^ Hz, z/syn,o,max = 4.0 x 10^"^ Hz, 
i^ssc,o,max = 1-0 X 10^^ Hz, F^y^.o = 3.5 X 10"^° erg cm'^ s'^, and Fss^.o = 4.5 x 10"^° erg cm'^ 
s~^. Substituting these values into the analytic estimate, we obtain 6 = 10.1, R = 3.5 x 10^^ 
cm, B = 0.55 G, 7br = 0.31 x 10^ 7max = 2.5 x 10^ Qe = 2.4 x lO'^ cm-^ s-\ and s = 1.8. 

In a similar way to the case of Mrk 421, we find that a satisfactory fit is obtained for 
the parameters as 5 = 11, i? = 1.0 x lO^^ cm, B = 0.20 G, 7i„ax = 2.0 x 10^ q^ = 1.7 x lO^^ 
cm~^ s~^, and s = 1.8. The thick solid curve in Figure |^ shows the predicted spectrum of 
Mrk 501 obtained by numerically calculating equations (|I|) and (^ for the best fit parameter 
set given above. The corresponding electron energy spectrum is shown by thick solid curves 
in Figure 0. The ratio of ub/uq derived from this numerical result is Me/mb = 22. 

As is for Mrk 421, the model predictions for the injection rates 50 times higher and 
lower than the best fit value are depicted by dashed and thin solid curves in Figures ^ and 
1^. The numerical values are tabulated in Table 2. Since the trend of uncertainties is the 
same as for Mrk 421, we do not repeat it here. 



4.1.3. PKS2155-304 

TeV emission from PKS 2155-304 {z = 0.117) was detected very recently by Durham 
Mark 6 Cherenkov telescope (Chadwick et al. 1999). Observables of PKS 2155-304 are 
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chosen as follows (Kataoka 2000); a = 0.2, z/syn,o,br = 4.0 x 10^^ Hz, z/syn,o,max = 2.0 x 10^^ 
Hz, i/ssc,o,max = 1-0 X lO^^ Hz, F^yn.o = 1-2 X 10"^ erg cm-2 s-\ and Fssc,o = 5.8 x 10~^° erg 
cm~^ s^^. Substituting these observables into the analytic estimate, we obtain 6 = 35.5, 
i? = 1.9 X 10^5 cm, B = 0.89 G, 7br = 0.34 x 10^ 7max = 0.76 x 10^ ge = 2.3 x 10"=^ cm'^ 
s~"^, and s = 1.4. 

In a similar way to the cases of Mrk 421 and Mrk 501, we find that a satisfactory 
fit is obtained for the parameter set given hj S = 33, R = 9.0 x 10^^ cm, B = 0.3 G, 
7max = 0.5 X 10^, ge = 2.8 X 10~^ cm~^ s~^, and s = 1.4. The thick solid curve in Figure ^ 
shows the predicted spectrum of PKS 2155-304 by numerically solving equations (|I|) and 
(0) for the best parameter set above. The corresponding electron energy spectrum is shown 
by the thick solid curve in Figure ^ The ratio of Ug/uB turns out to be Ug/uB = 3 for this 
model parameters. 

As is for Mrk 421 and Mrk 501, the model predictions for injection rates 50 times higher 
and lower than the best fit value are depicted by dashed and thin solid curves in Figures || 
and ^ The numerical values are tabulated in Table |^, and the trend of uncertainties is the 
same as for Mrk 421 and Mrk 501. 



5. SUMMARY AND DISCUSSION 

To summarize, within the framework of the one-zone synchrotron self-Compton model, 
we determined the numerical values of the physical quantities of TeV blazars Mrk 421, Mrk 
501, and PKS 2155-304 in quiescent states. Those values are searched for by solving the 
kinetic equations of electrons and photons taking proper account of injection, escape, and 
cooling of electrons and by comparing predicted radiation spectra with observations. The 
best fitted parameters and uncertainties are estimated. It is shown that the ratio of the 
energy density of electrons to that of magnetic fields can be determined within a factor 
of a few and the ratio is about 5, 22, and 3 for Mrk 421, Mrk 501, and PKS 2155-304, 
respectively. Thus, the emission region of TeV blazars is particle dominated. For other 
parameters, 6 and 7max are also determined within a factor of a few, while the magnetic field 
strength and the size have an uncertainty of a factor of 10 and 30, respectively. Since these 
results have important implications for the fundamental understanding of the formation 
and bulk acceleration of relativistic jets, below we discuss some of the further issues to be 
explored. 
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5.1. Ue/uB Ratio 

First, we discuss the value of 7min because it is important for determining the Ue/u-Q 
ratio. From the theoretical standpoint, in the case of pair plasma jets, we regard it most 
likely that 7min ~ T, since the shock first thermalizes a bulk population of particles and 
then accelerates them from this pool. On the other hand, for the normal electron-proton 
plasma, there is a wide range of possibilities about 7min- One extreme case is that protons 
and electrons are separately thermalized, which leads to 7min ~ T- Such a separate 
thermalization is supposed to be realized for non-relativistic shocks in supernova remnants. 
The other extreme case is that electrons and protons attain an equilibrium state, which 
means 7min ~ V^y . This picture is conventionally assumed for models of cosmic 7-ray 
bursts (Sari, Piran, & Narayan 1998). In this case 7min ~ 10^ may be realized, although the 
real value is likely to be between these two values. 

From the observational standpoint, it is also difficult to determine 7min- Since the 
observed synchrotron frequency from an electron with the Lorentz factor 7min is given by 

^sy„,o,min = 1-2 X 106^7^.^-^ , (34) 

electrons with 7jnin ~ 10 and 10'^ emit synchrotron photons of ~ 10® and 10^^ Hz, 
respectively, for typical values B ~ O.IG and 5 ~ 10. The former is well below the 
self-absorption frequency and the emission in the low frequency band is dominated by that 
from more extended regions. Since the latter frequency is well above the self-absorption, 
we must modify the model such that emission below lO^'^ Hz is not from X-ray emitting 
one-zone region but from a separate region. 

In our one-zone model, to attain the conventional equi-partition state of Ue = ub by 
simply changing 7inm alone, we need to adopt 7inin ~ 1.7 x 10^, 3.2 x 10"^, and 0.5 x 10^ 
for Mrk 421, Mrk 501, and PKS 2155-304, respectively. It is of some interest that these 
analytically estimated value is near to 10^ mentioned above. For 7min ~ 10^, the predicted 
spectra around near infrared and optical bands are difficult to match the observed data 
in the one-zone model. In this case, whole synchrotron spectra may be reproduced by a 
superposition of emission from inhomogeneous jets. An example of this kind of modeling is 
shown in the recent work of [Katarzyhski, Sol, fc Kus(2001)| about the broad-band spectra 
of Mrk 501 in the flaring stage by a one-zone SSC blob in a conical jet. The blob explains 
X-ray and 7-ray bands, while the jet explains the spectrum from optical to radio band. We 
note that their Ue/ub ratio for the blob is similar to ours, irrespective of the actual value of 
7min) for this flaring stage of Mrk 501 although they did not explicitly mentioned. 
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5.1.1. The Case of •yrnin = 10^ 

As discussed above, it is important to examine if the alternative case of 7min ~ lO'^ and 
Ue/u-Q <^ 1 for normal plasma can reproduce the high energy part of the emission spectra in 
quiescent states. 

First, we check the case where 7br is still larger than 7min- The other parameters are 
fixed as the best fitting ones obtained in the previous section. The numerical results for 



Mrk 421 are shown in Figures [1^ and |TT]. As is seen in Fig. |T0|, the low energy part of 
both synchrotron and inverse Compton emission is short of the observation. In this case, 
the numerical value of Ue/u-Q turns out to be 2 and yet the emission region is kinetic power 
dominated. Hence we can rule out the possibility of Uc/ub < 1 for 7min < 10 



4 



The second case is 7min > 7br for which we analytically examine the expected trend. 
When 7br < 7min is satisfied, the electron energy spectrum is given by 

^^6(7) = ge4,esc7mfn^Sbr7~^ for 7br < 7 < 7min 

ne(7) = ge4,esc7br7~''~^ for 7min < 7 < 7max- (35) 



This regime is called fast cooling by |Sari, Piran, fc JNarayan(1998)| in the models of gamma 



ray bursts. In this case, instead of Eq. ( PE| ) we obtain 

^_Q -1 -1 f, , /ssc ^ /ssc 7br ln(7min/7br) /ocN 

«B - ' '^ '^ V^ fsyn ) /s.n 7^n^ C: ^^^'d^ ' ^ ^ 

Since the combinations of the observables z/syn,o,max, '^ssc.o.max, -Z^syn.o, and Lssc,o are the same 
as before, from Eqs. (^, (|^, and (P^T]), 57niax, -B(^7max5 ^^^ 6'^R^B^ are not changed. From 
Eq. (|T2p, 7br is proportional to R^^B~^ if we neglect Compton cooling. Thus, we obtain 
B (X 6 (X 7br, R oc 7^1.'^, and 7max oc %^^. From this, we see that stronger cooling, i.e., smaller 
values of 7br, implies smaller magnetic field strength and beaming factor and larger size 
and 7max- This is not a favorable choice of parameters, because their direction of changes 
is contradictory to the observational facts about strong beaming and rapid time variability. 
Therefore fast cooling regime is not favored. We thus conclude that our conclusion on 
kinetic power doninance is fairly robust. 



5.2. Proton Components 

It is interesting to know what constraints on the jet material can be obtained from this 
analysis. Let us assume that the jet consists of relativistic electrons and cold protons and 
that there are neither thermal electrons nor relativistic protons. Then, using the electron 
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number density for 7jnin = 10 and charge neutrality, the kinetic power of cold protons is 
estimated to be 2.4 x 10^^ erg s'^ 7.8 x 10"^^ erg s'^ and 4.0 x 10"^^ erg s'^ for Mrk 421, 
Mrk 501, and PKS 2155-304, respectively. These values are several to 10 times larger than 
the kinetic power of electrons. This small ratio of the kinetic power of protons to that of 
electrons, not withstanding the large mass ratio, is due to a large average Lorentz factor 
of relativistic electrons. Because the proton kinetic power does not exceed the Eddington 
luminosity for a representative black hole mass in AGN, we do not make strong argument 
about jet material in this way alone. Also, in the case for 7min = 10^, proton kinetic power 
is less than that of relativistic electrons. Estimation of the large scale kinetic power of these 
sources can in principle discriminate between these two possibilities, i.e., proton-electron 
jets or electron-positron jets. Considering the weakness of the extended radio emission of 
TeV blazars, the same kind of analysis for GeV blazars seems to be more promising. We 
will analyze GeV blazars in the future. 



5.3. Further Comments 

We discuss the injection index s. For TeV blazars studied in this work, s is smaller 
than 2. Some recent work have reported that in Fermi acceleration at ultra-relativistic 
shocks s is larger than 2 ([Gallant, Achterberg, and Kirk(1999) ; Kirk et al.(2000)| ), different 



from our adopted values based on observed spectra. This is an interesting open question 
and future research is awaited. 

As for the ratio of the energy densities of the relativistic electrons to the magnetic 
fields, what is physically more meaningful may be the ratio (^e + u^a.d/C3)/uB, because this 
corresponds to the ratio of the injected kinetic power to the Foynting power, while Us/ub 
corresponds to the remaining ratio after radiative cooling. The former ratio becomes 10, 41, 
and 6 for Mrk 421, Mrk 501, and PKS 2155-304, respectively. Thus, the conclusion does 
not change much, although somewhat strengthened. 

As most of the previous work, we also neglected the correction for the absorption of 
TeV 7-rays due to GIB. For Mrk 421 and Mrk 501 in quiescent states, this neglect seems 
safe, if the observed spectra do not much extend over the TeV range. However, for PKS 
2155-304 of the redshift of 0.117, the optical depth for 0.35 TeV photons amounts to 
about 0.5, although the exact value depends on the GIB intensity and spectrum. Thus the 
intrinsic SSG luminosity of PKS2155-304 should be a little higher than the estimate given 
in this paper. Consequently, the ratio Uq/ub is somewhat higher than the value obtained 
above and our conclusion of the kinetic power dominance is strengthened. 
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Now, we shortly compare our results with previous work of others, [lavecchio et al. 



1998)1 derived values of physical parameters in a similar way to our analytical method. 
However, in their paper, electron energy distribution is not calculated self-consistently. 
It means that normalization, index, and 7br of relativistic electrons are not derived from 
solving the cooling and injection processes but obtained only from a fitting of the observed 
spectrum with double power-law spectrum. Here we examine the case of Mrk 421, because 
this is the only source for which multi-frequency spectral fitting was done and the number 
density of the accelerated electrons was estimated by them. 

Their best fit parameters for Mrk 421 are B = 0.15 G, 6 = 25, and R = 2.7 x 10^^ 
cm, and their electron energy spectrum is given by ne(7) = K'y~"'^{l + 7/7br)"'^~"'^5 with 
rii = 2.2, 77-2 = 4.5, 7br = 5.6 x 10^, and K = 1.7 x 10^. Note that they are significantly 
different from our values. Their electron spectral index 2.2 is steeper than 1.6 in this work, 
and it produces a much steeper synchrotron spectrum at low energies. Their estimate of 
the emission region is more compact and more strongly beamed than ours. Their parameter 
values correspond to Ue/uB = 371 for a choice of 7min = 10 and 7niax = 10^, and the ratio 
does not much change for a different choice of 7min and 7max- Thus, the adopted values by 
Tavecchio et al. (1998)| would mean even more kinetic power dominated states than our 



results. Major reason for this difference is due to the choice of s, because replacing s = 2.2 
by s = 1.6 leads to Us/ub = 19. 

Although we do not make further comparison with other work, we again emphasize 
that our method is superior to previous ones in that the break energy of electrons is 
self-consistently determined and direct spectral fitting is made. 

We thank an anonymous referee for useful comments which help us to improve this 
paper. We are grateful to J. Kataoka for kindly providing us with observational data and 
for useful discussions. This work is supported in part by a Grand-in-Aid from Ministry of 
Education, Science, Sports and Culture of Japan (11640236, F. T.). 
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Table 1. Mrk 421 Physical Parameters from SSC Analysis 







high injection 


best fit 


low injection 


5 




15.9 


12 


9.1 


R (cm) 




5.3 X 10^^ 


2.8 X 10^6 


1.4 X 10^^ 


5(G) 




0.44 


0.12 


0.036 


O'max 




8.3 X 10^ 


1.5 X 10^ 


3.45 X 10^ 


S 




1.6 


1.6 


1.6 


Qe (cm"2 S"^)"" 




4.9 X 10-^ 


9.6 X 10"^ 


1.9 X 10-^ 


nf* (cm-3 Y'^ 




1.1 X 10^ 


1.1 X 10^ 


1.1 


b)" 




2.3 X 10^ 


3.1 X 10^ 


4.3 X 10^ 


'77 




8.9 X 10-3 


6.9 X 10-3 


5.5 X 10-3 


^yn,o (erg s"^) 




1.3 X 10^^ 


1.3 X 10^5 


1.4 X 10^5 


Lssco (erg s"^) 




0.5 X 10^5 


0.7 X 10^5 


0.9 X 10^5 


-^syn+ssc,o (,erg S 


-') 


1.8 X 10^5 


2.0 X 10^5 


2.3 X 10^5 


Lpoy (erg s-i) 




6.9 X 10^2 


8.1 X 10^2 


1.0 X 10^3 


^ckin (erg s"^) 




1.8 X 10^3 


4.0 X 10^3 


7.8 X 10^3 


-t^syn+ssc (erg S 


') 


7.2 X 10^2 


1.4 X 10^3 


2.8 X 10^3 


^p.kin (erg s~^)'^ 




1.5 X 10^^ 


2.4 X 10^^ 


3.3 X 10^^ 


-^Cjkiii/ -i^poy "^c 


'Jub 


3 


5 


8 



^7min is fixed at 10. 



b„tot 



and (7) are the total number density and the average Lorentz factor 
of relativistic electrons, respectively. 

^uLi, ^ 1.3 X 10^^ erg s~^ at 0.3 keV is adopted in calculating r^^. 

'^-^p.kin is calculated assuming that the cold proton number density is the 
same as that of relativistic electrons. 



-21- 



Table 2. Mrk 501 Physical Parameters from SSC Analysis 





high injection 


best fit 


low injection 


5 


14.5 


11 


8.3 


R (cm) 


2.0 X 10^^ 


1.0 X 10^6 


5.3 X 10^6 


5(G) 


0.69 


0.20 


0.059 


O'max 


9.0 X 10^ 


2.0 X 10^ 


3.8 X 10^ 


S 


1.8 


1.8 


1.8 


Qe (cm"2 S"^)"" 


8.7 X 10-2 


1.7 X 10-3 


3.4 X 10-5 


n*°* (cm-3 )'^'b 


3.4 X 10^ 


3.4 X 10^ 


3.6 X 10^ 


(7)" 


1.0 X 10^ 


1.2 X 10^ 


1.5 X 102 


'77 


2.2 X 10-2 


1.8 X 10-2 


1.2 X 10-2 


^yn,o (erg s"^) 


8.3 X 10^^ 


7.8 X 10^^ 


8.2 X 10^^ 


Lssco (erg s"^) 


9.9 X 10^^ 


1.0 X 10^5 


1.2 X 10^5 


-t^syn+ssc,o (erg S ) 


1.8 X 10^5 


1.8 X 10^5 


2.0 X 10^5 


Lpoy (erg s-1) 


2.0 X 10^2 


2.4 X 10^2 


3.4 X 10^2 


^ckin (erg s"^) 


3.0 X 10^3 


5.2 X 10^=^ 


1.1 X 10^^ 


-t^syn+ssc (erg S ) 


8.6 X 10^2 


1.5 X 10^=^ 


2.9 X 10^3 


^p.kin (erg s~^)'^ 


5.2 X 10^^ 


7.8 X 10^^ 


1.3 X 10^5 


J^c,km/ -l^poy = '^c/'^B 


15 


22 


32 



^7min is fixed at 10. 



b„tot 



and (7) are the total number density and the average Lorentz factor 
of relativistic electrons, respectively. 

^uLj, ^ 7.9 X 10^3 erg s-^ at 0.3 keV is adopted in calculating t^^. 

'^-^p.kin is calculated assuming that the cold proton number density is the 
same as that of relativistic electrons. 
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Table 3. PKS 2155-304 Physical Parameters from SSC Analysis 







high injection 


best fit 


low injection 


5 




43.6 


33 


25.0 


R (cm) 




1.9 X 10^^ 


9.0 X 10^5 


4.5 X 10^6 


5(G) 




1.1 


0.30 


0.089 


O'max 




2.3 X 10^ 


5.0 X 10^ 


9.8 X 10^ 


S 




1.4 


1.4 


1.4 


Qe (cm"2 S"^)"" 




1.4 X 10-3 


2.8 X 10-5 


5.6 X 10-^ 


nf* (cm-3 Y'^ 




2.5 X 10^ 


2.4 X 10^ 


2.3 


b)" 




2.9 X 10^ 


4.8 X 10^ 


7.3 X 10^ 


'77 




5.4 X 10-3 


4.6 X 10-3 


3.7 X 10-3 


^yn,o (erg s"^) 




3.4 X 10^6 


3.2 X 10^6 


3.4 X 10^6 


Lssco (erg s"^) 




1.1 X 10^6 


1.5 X 10^6 


2.0 X 10^6 


-^syn+ssc,o (,erg S 


-^) 


4.5 X 10^6 


4.7 X 10^6 


5.4 X 10^"^ 


Lpoy (erg s-i) 




4.1 X 10^3 


4.0 X 10^3 


4.9 X 10^3 


^ckin (erg s"^) 




5.2 X 10^3 


1.0 X 10^^ 


2.2 X 10^^ 


-t^syn+ssc (erg S 


') 


2.4 X 10^3 


4.3 X 10^3 


8.7 X 10^3 


^p.kin (erg s~^)'^ 




3.3 X 10^^ 


4.0 X 10^^ 


5.7 X 10^^ 


-^Cjkiii/ -i^poy "^c 


'Jub 


1 


3 


4 



'^Tmin is fixed at 10. 

b^tot ^j^^ ^^^ g^j,g ^]-^g total number density and the average Lorentz factor 
of relativistic electrons, respectively. 

^uLi, ^ 4.1 X 10^5 erg s~^ at 0.3 keV is adopted in calculating r^^. 

'^-^p.kin is calculated assuming that the cold proton number density is the 
same as that of relativistic electrons. 
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Fig. 1. — Schematic picture of the multi-frequency spectrum of a typical TeV blazar. 
Here Lsyn,o and Lssc,o are the observed bolometric luminosities of synchrotron and SSC 
components, respectively. Corresponding to the break in the relativistic electron energy 
spectrum in the emission region, a break feature appears in the observed synchrotron 
spectrum. Around the TeV energy region, the Klein-Nishina effect suppresses the observed 
flux compared to the Thomson regime. We exclude the information of z/ssc,o,br, because of 
sparse observational data points and complication from the Klein-Nishina effect. 
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Fig. 2. — Schematic picture of the relativistic electron energy spectrum. At high Lorentz 
factors, radiative coohng decreases the number density of electrons and leads to a break in 
the spectrum. 
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Fig. 3. — Ratio u^/u-q in TeV blazars. This is calculated according to equation (pHI). For 
a TeV blazar u^ is typically about one order of magnitude larger than u-q. We confirm this 
analytic result by numerical calculations. The discrepancies between this analytic estimation 
and numerical results are at most within a factor of a few. 
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Fig. 4. — One- zone SSC model spectra for the steady state emission of Mrk 421. The thick 
solid line shows the best fit spectrum, where adopted parameters are 5 = 12, R = 2.8 x 10^^ 



cm, 5 = 0.12 G, 7„ 



1.5 X 10^, q^ 



9.6 X 10-6 cm-3 s-\ 



s = 1.6, and u^/ub = 5. The 



dotted line shows the spectrum obtained using the analytic estimates for Mrk 421. The thin 
solid and dashed lines show the spectra of low and high injection models, respectively, to 
indicate the uncertainty range of the spectral fitting. 
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Fig. 5. — Electron energy spectrum and kinetic power of Mrk 421 corresponding to Figure 
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Fig. 6. — One- zone SSC model spectra for the steady state emission of Mrk 501. The thick 
solid line shows the best fit spectrum where adopted parameters are 5 = 11,/?= 1.0 x 10^^ 



cm, B = 0.20 G, 7^ 



2.0 X 10^ gc 



1.7 X 10 ^ cm ^ s" 



s = 1.8, and u^/u-q = 22. The 



dotted line shows the spectrum obtained using the analytic estimates for Mrk 501. The thin 
solid and dashed lines show the spectra of low and high injection models, respectively, to 
indicate the uncertainty range of the spectral fitting. 
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Fig. 7. — Electron energy spectrum and kinetic power of Mrk 501 corresponding to Figure 
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Fig. 8. — One-zone SSC model spectra for the steady state emission of PKS 2155-304. 
The thick solid line shows the best fit spectrum where adopted parameters are 5 = 33, 
i? = 9.0 X 10^5 cm, B = 0.30 G, 7max = 0.5 x 10^ q^ = 2.8 x lO'^ cm^^ g-i^ ^ ^ ^4^ ^^^^ 
Uc/ub = 3. The dotted line shows the spectrum obtained using the analytic estimates for 
PKS 2155-304. The thin solid and dashed lines show the spectra of low and high injection 
models, respectively, to indicate the uncertainty range of the spectral fitting. 
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Fig. 9. — Electron energy spectrum and kinetic power of PKS 2155-304 corresponding to 
Figure ^. 
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Fig. 10. — One-zone SSC model spectra for the steady state emission of Mrk 421 for 
7min = 1 X 10^, and Uq/ub = 2. The thick solid line shows the spectrum for this case; 
other parameters are the same as Figure ^. The thin solid line shows the best fit spectrum 
of Mrk 421 shown in Figure H. 
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Fig. 11. — Electron energy spectrum and kinetic power of Mrk 421 corresponding to Figure 



